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The  Calculation  of  the  Heat  Transfer  Associated  with  the 
Injection  of  Methane  Gas  at  the  Stagnation  Point  of  the 
Dissociated  Boundary  Layer  of  the  Laminar  Air  Flow 

Yan  Longhua  Yu  Xintian 

(Beijing  Institute  of  Aerodynamics) 

ABSTRACT 

By  using  methane  gas  as  the  injecting  medium,  this  paper 
examines  in  detail  the  effect  of  the  chemical  combustion 
reaction  on  the  heat  transfer  at  the  stagnation  point  of  the 
laminar  flow  boundary  layer  associated  with  a  blunt  rotating 
object.  The  following  assumptions  are  made  concerning  the 
characteristic  transport  properties  of  the  mixed  gases.  The 
Pr  and  Le  numbers  as  well  as  are  set  to  be  constants. 

The  boundary  layer  is  assumed  to  be  in  a  state  of  thermodynamic 
equilibrium.  Some  calculations  are  also  conducted  by  treating 
the  methane  gas  as  an  inert  gas  so  that  the  effect  of  heat 
release  associated  with  the  chemical  combustion  reaction  can 
be  examined  by  cross  comparison. 

I.  INTRODUCTION 

After  the  appearance  of  the  thermal  ablation  phenomenon 
at  the  front  section  of  the  material,  the  thermal  environment 
around  the  front  section  of  a  supersonic  flying  object  changes 
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drastically.  The  change  occurs  since  the  speed  in  which 
foreign  matter  is  introduced  into  the  boundary  layer  is  high. 

This  can  decrease  the  heat  flow  added  on  the  surface  of  the 
object  by  its  environment.  In  the  same  time,  however,  the 
heat  flow  can  also  be  increased  by  this  added  matter  if  it 
undergoes  chemical  reaction  which  releases  heat  within  the 
boundary  layer.  Extensive  research  has  been  conducted  abroad 
concerning  the  computation  of  heat  flow  associated  with  the 
chemical  reactions  of  added  matters  in  the  laminar  flow  boundary 
layer.  For  example,  references  [ll  -  [3~)  addressed  this  problem 
by  solving  the  solution  numerically.  They  studied  the  thermal 
ablation  of  graphite  by  solving  the  equations  of  the  stagnation 
point  of  the  laminar  flow  boundary  layer  in  which  multiple 
sets  of  elements  and  combustion  reactions  exist.  Based  on 
their  results,  it  seems  that  even  though  there  are  many 
chemical  reaction  equations  which  increase  the  number  of 
the  system  elements  and  the  transport  coefficients  of  the  mixed 
gases  are  also  increasingly  complicated,  theoretical  study 
can  still  be  done  with  some  reasonable  or  some  fixed  simplifying 
assumptions.  The  theoretical  study  can  provide  insight  into 
the  mechanism  of  the  chemical  reactions  which  affect  the  thermal 
environment.  Reference  [4}  assumed  that  the  transport  coefficients 
9  1,  Pr  =  l  and  Le  =  l.  It  applied  the  combustion  area 
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model  and  derived  numerical  solutions  to  examine  separately 
the  combustion  processes  within  the  boundary  layers  of  a 
blunt  half  sphere  and  a  flat  plate,  with  hydrogen  gas  as 
the  injecting  medium.  In  the  area  of  the  approximated  solutions, 
reference  [5)  proceeded  with  some  general  discussion  and  provided 
some  qualitative  concepts  about  the  equilibrium  and  the  frozen 
catalytic  wall  associated  with  the  chemical  reactions  of  the 
added  matter  at  the  stagnation  point  of  the  laminar  flow 
boundary  layer  by  assuming  Le=l.  Similar  discussion  was  also 
carried  out  under  the  frozen  catalytic  wall  condition  of  Le  1. 
For  certain  materials,  such  as  the  mixed  gases  with  Le  numbers 
close  to  1,  it  also  provided  the  methods  for  obtaining 
approximated  quantitative  estimates.  Reference  £6}  treated  the 
carbon  combustion  problem  by  using  the  combustion  surface  model. 
This  model,  however,  lajjks  generality. 

The  necessary  results  obtained  from  numerical  computations 
are  still  incomplete  as  far  as  the  research  of  the  engineering 
computational  method  is  concerned.  This  is  why  that  this  paper 
considers  the  use  of  methane  gas  as  the  injecting  medium  for 
the  numerical  computation  of  the  combustion  of  the  added  matter 
at  the  stagnation  point  of  the  boundary  layer.  The  computation 
assumes  that  after  methane  is  injected  into  the  high  temperature 
boundary  layer,  the  following  two  major  chemical  reaction 
processes  and  a  decomposition-reconstitution  reaction  will  occur* 
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CH.  +  3Q./2  CO  +  2H,G 

(1) 

CO  +  CV2  —  CO, 

(2) 

O,  — •  20 

(3) 

* 

This  paper  was  presented  at  the  Huang  Mountain  Thermal 

Conducting  Material  Science  Meeting  in  October,  I98I. 

When  the  static  temperature  is  not  much  higher  than  4000K, 
we  can  reluctantly  neglect  the  effect  of  the  decomposition 
of  the  nitrogen  molecule.  H20  will  not  be  decomposed  into 
hygrogen  and  oxygen  and  all  other  possible  chemical  reactions 
will  be  treated  as  of  only  secondary  importance  and  hence 
neglected.  At  the  stagnation  point,  the  velocity  is  zero  so 
we  can  assume  that  the  chemical  reaction  mentioned  above  will 
be  at  equilibrium  state.  The  method  for  computing  the  transport 
coefficient  of  some  element  (such  as  methane)  is  unknown  so 
we  assume  that  1,  Pr  and  Le  to  be  constants.  The  limited 
computational  examples  use  the  charac"  eristic  conditions  of 
a  flow  temperature  T_  =  5000 K  at  the  stag"  Ion  point  and  a 
pressure  P  which  equals  one  atmosphere. 

II.  BASIC  EQUATIONS  AND  TREATMENT  METHOD 

Starting  from  the  equations  of  continuity,  element 
diffusion  and  energy  (13)-(15)  as  well  as  the  equation  of 
motion  (1?)  associated  with  the  basic  equations  of  a  laminar 
flow  boundary  layer  of  a  rotating  body  in  decomposed  gases 
derived  by  reference  [7]  ,  we  can  apply  the  following 
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(4) 


transformation  by  neglecting  the  thermal  diffusion 

(  —  \ 

and  ws  can  introduce  the  dimensionless  variables  ft-u/u0. 
g=  (htu72)hg,  f  =  l  f?  +  f(  0)  f  0=T/Tg.  Using  a  method  similar 
to  that  in  reference  [7]  *  we  can  further  use  the  conditions 
of  the  air  flow  at  the  stagnation  point  and  the  shape  of  the 

r-  x,  «;«*,,  ut==  (du./<tr)j  ; 

object/to  convert  the  equations  mentioned  above  into  a  set 
of  ordinary  differential  equations  with  dimensionless  coefficients 
(this  will  automatically  be  satisfied  after  the  transformation 
of  the  continuity  equation) t 

l(/Le/PrK„),  +  /C„  - 

#“1,2, — 7  (5) 

K//Pr)*,J,  +  fg ,  +  <(/<Le  -  l)/Pfi  HSU  -  *?)C„]},  -  0  (6) 

(U„),  +  //„  +  (*>-«)/ 2-0  (  ?  ) 

In  these  equations  the  subscripts  "  £ "  represent  the  derivatives 
taken  with  respect  to  £  •  Here  we  have  already  assumed  that 
the  diffusion  coefficients  of  all  the  system  elements  are 
the  same.  In  other  words,  the  Lewis  numbers  are  the  same.  Other 
symbols  include  the  enthalpy  h “Zc^(h^-h^°) ,  h^=JT  cp^  dT, 

COi  -  the  mass  formation  rate  of  the  i-th  system  element  in 
the  chemical  reaction,  T  -  temperature,  density  and 

-  the  mass  concentration  of  the  i-th  system  element.  The 
chemical  symbols,  molecular  weights  and  the  formation  enthalpy 
h^  (also  called  the  reaction  heat)  of  the  i  system  elements 
are  shown  in  Table  1.  hg  is  the  total  enthalpy.  The  subscript 
"c"  indicates  the  condition  at  the  outer  fringe  of  the  boundary 
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layer  while  the  subscript  "s"  indicates  the  condition  at  the 
stagnation  point.  The  meaning  of  the  two  are  the  same  here. 
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Table  1.  The  molecular  weights  and  the  formation  enthalpy 
(h?  ,  cal/kg)  of  the  constituent  elements  within  the  boundary 
layer. 

1.  constituent  element  i  2.  chemical  symbol 

3.  molecular  weight 

We  will  now  further  assume  that  Le  (-^>Dcp/^  )  is  a 
constant.  This  is  different  from  the  actual  condition.  The 
differences  between  the  maximum  values  and  the  minimum  values 
of  the  diffusion  coefficients  associated  with  different 
constituent  elements  can  be  quite  large,  sometimes  exceeding 
b0%.  The  precise  computation  has  two  difficulties!  the 
computational  method  for  the  calculation  of  our  desired  gas 
transport  coefficients  at  high  terperatures  is  still  being 
studied  and  the  equations  for  the  computation  of  the  transport 
coefficients  are  extremely  complicated.  The  solutions  of  the  main 
equation  sets  will  be  difficult  to  obtain  if  we  were  to  consider 
this  problem.  It  will  take  several  more  computations 


and  it  is  not  advantageous  for  any  preliminary  study.  Similarly, 


we  assume  that  =  constant,  with  ^  being  the  viscosity 

coefficient.  We  also  assume  that  Pr  =•  constant,  with  Pr  t 

X.  is  the  thermal  conductivity  coefficient  ,  "c  -  Tc.  c  . 

p  1  pi 

and  Jj =  In  order  to  compensate  for  the  deficiency 
of  the  |s  1  assumption,  we  have  selected  several  different 
constants  for  Q  in  the  computation.  We  have  also  treated  Pr 
and  Le  in  a  similar  fashion. 

We  assume  that  there  are  three  chemical  reactions  occurring 
simultaneously  within  the  boundary  layer  and  they  all  reach 
chemical  equilibria.  There  are  a  total  of  seven  constituent 
elements  and  the  nitrogen  gas  does  not  participate  in  the  chemical 
reaction.  In  this  way,  each  constituent  element  should  satisfy 
three  mass  reaction  principle  relationship  equations  corresponding 
to  the  three  chemical  reactions  ( 1 ) — ( 3 )  as: 


C?CJ  - 

Kn  U,M\ 

(8) 

C,C\* 

uy  mmX' 

C,  _ 

(  Q\ 

cZF 

MMV 

\  7  / 

£»- 

KrU\ 

(10) 

C, 

Here  Mp  =  2^  •  These  three  relationship  equations 

can  substitute  the  three  equations  in  (5). 

Within  the  boundary  layer  which  contains  the  seven 
constituent  elements,  as  long  as  the  mass  flow  is  fixed, the 
four  chemical  elements  H,C,0  and  N  will  be  conserved  regardless 
of  the  condition  of  the  chemical  reaction.  This  is  due  to  the 
principle  of  conservation  of  mass  for  a  closed  system.  We  will 
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now  introduce  four  new  mass  quotients  of  similar  chemical 
elements  as  j  'C"1  =;  4-  =  Gc  =  C2 

clM2/Mi  +  C3M2/M3j  ^  =  'cQ  W  C5  4-  2C2M5/M2  1.5  4C6; 

Cy  =  Cjy=  C^,.  The  mass  formation  rate  corresponding  to  these 
four  values  are  zero.  We  can  thus  obtain  fov  relationship 
equations  as.-a)^  4-  0.5  60^  =0;  ^2 +^1M2/  +  —  °* 

60^+  2  Cv)2M3/M2  H-1.5  +CJg  =  0;  (aJ^  .We  can  now 

combine  the  diffusion  equation  (5)  according  to  the  fashion 
of  these  four  relationship  equations  and  obtain  four  more 
new  equations. 

[(/WPr)C„J,  +  fCn  -  0 

,  -1,2,5, 7  (11) 

The  general  solutions  of  these  four  equations  are  similar. 

We  can  establish  three  relationship  equations  as 

ft  ~  ft~  -  -  ft—  £»  ~  -  <**)  (12) 

The  subscript  MwM  here  indicates  the  condition  at  the  surface  of 
the  wall.  By  differentiating  these  equations  with  respect  to  £ 
and  applying  the  relationship  between  the  corresponding  boundary 
conditions,  we  can  transform  them  into  the  following  three 
new  equations. 

CJU ,  +  CJ «,  -  C./(2*#.) 


C,  +  C,  +  c,  +  C.  +  C,  +  C,  4-  Ct  —  1 

(13) 

£3.  +  Si-  +  S±.  +  +  ^7 

(14) 

M,  2  M.  2Af,  W,  21#. 

CfcMj  C*  _  0.2095  C- 

07108 

(15) 

We  can  use  these  three  equations  to  substitute  any  three  equations 
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in  (11).  We  will  select  and  retain  the  equation 

l(fLe/Pr)C„],  ■+■  fCm  ™  •  (16) 

We  can  ultimately  obtain  the  numerically  solvable  ordinary 

differential  equations  (6),  (7),  (16)  and  the  algebraic 

equations  (13)-(15)»  (8)-(10).  After  going  through  the 

transformation  (4),  we  can  write  the  corresponding  boundary 

conditions  of  the  differential  equations  as  : 

1-0:/,  -  0, 

/ -  (A) 

/  -/.  — un.V.’M.w*,'*)  »  - 

/.(  W,  V,  -  C„)  -  (/Le/PrXC*). 
q  — *  x  :  /,  — ►  I »  g  0 

Under  an  outer  flow  temperature  of  5000K,  we  should 
consider  the  contributions  of  the  molecular  translation, 
rotation  and  degrees  of  freedom  of  oscillation  on  the  isobaric 
specific  heat  of  the  i  constituent  elements.  The  effect  of  the 
electron  excitation  energy  can  be  neglected.  In  order  to 
simplify  the  computational  equations,  we  will  use  the  following 
approximated  equations  with  errors  less  than  Z/o  : 


P^ 


c „  -  2.4894  -  0.60869(0-"  ’  +  0.50) 
t„  -  0.41032  -  0.04936(0-*’  +  0.50) 
c„  -  0.37966  -  0.04 166(J“*-’  +  0.50) 

-  1.003  -  0  051310  -  0  202630-*-’  (C) 

-  0.31319  -  0.02319(0-°  *  +  0.50) 

-  0.31 1 

e„  -  0.3856  -  0.045846(0-*-’  +  0.50) 


In  this  way  we  can  compute  c^  .  From  the  three  dimensionless 
parameters  obtained  from  the  combination  of  the  transport 
coefficients  of  the  mixed  gases  ,  we  now  select  certain  combinations 
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of  the  following  values:  set  ^  to  be  0.5* 1,1. 5  and  2;  set 
Pr  to  be  0.5*  0.71  and  1;  set  Le  to  be  1,  1-4,  2  and  4. 

Based  on  the  assumptions  mentioned  earlier,  we  will  take 
l  =  1,  Pr  =  0.71  and  Le  —  1.4  as  the  basic  condition. 

We  then  apply  the  computational  method  of  thermal 
chemistry  to  obtain  the  chemical  equilibrium  constants 
associated  with  the  chemical  reaction  equations  (l)-(3). 

In  order  to  simplify  the  computations,  we  will  use  the  following 
approximated  equations  at  temperatures  below  8000K  i 

-  27*35 +  3.9056 

%  K, ,  -  1 34*3  (  TjB)  -  4.39 1 9  <  D  * 

la*:,,  -  16.354  -  61530/(T,fl) 


We  can  now  numerically  compute  the  ordinary  differential 
equations  along  with  their  boundary  conditions  and  the  algebraic 
equations. 

We  can  obtain  the  heat  flow  transferred  into  the  object 
based  on  the  energy  equilibrium  between  the  upper  and  the  lower 
surfaces  of  the  object  K).  After  going  through  transformation  (4), 
we  can  then  obtain  a  heat  flow  computational  equation  similar 
to  that  in  reference  C7}  • 


The  relationship  equation  for  the  frictional  coefficient  is: 

c,KV  2  -  \/T//„(0),  R„  -  *■*/*..  (F) 


III.  DISCUSSION  OF  THE  RESULTS  OF  THE  COMPUTATIONS 


10 


A 


The  conditions  of  having  a  temperature  of  5000K  at  the 
outer  fringe  of  the  stagnation  point  boundary  layer  and  a 
stagnation  point  pressure  of  1  atmosphere  are  comparable 
to  the  conditions  of  an  object  flying  at  an  altitude  of 
thirty  thousand  meters  with  a  Mach  number  of  15.  We  will 
now  assign  exclusively  the  "basic  condition"  as  the  variation 
in  the  properties  of  the  stagnation  point  boundary  layer  when 
the  dimensionless  parameter  -fw  increases  from  0  to  0.533 
under  the  conditions  that  i=l,  Pr  =  0.?l,  Le  =  1.4  and  0(0)  =0.4. 
We  will  now  proceed  with  the  discussion  based  on  five  different 
types  of  local  modifications  of  the  "basic  condition". 

1.  Basic  Condition  In  the  absence  of  the  injection, 

the  heat  flow  at  the  stagnation  point  of  the  laminar  flow 
boundary  layer  can  be  computed  based  on  equation  (63)  in 
reference  [7]  .  We  can  obtain  qQ  =  1816  CfcMucx)s3  • 

With  a  small  amount  of  methane  leaking  from  inside  the  wall, 
the  methane  will  bum  on  the  surface  of  the  wall  instantly 
since  the  wall  temperature  is  higher  than  the  ignition 
temperature  of  mpthane.  The  heat  released  from  this  process 
will  be  added  directly  on  the  surface  of  the  wall,  causing 
the  transfer  of  heat  flow  to  the  surface  of  the  wall  to  increase 
sharply  with  the  increase  in  -fw.  This  kind  of  phenomenon  can 
be  observed  in  Figure  1.  For  example,  when  -fw=0.02, 
q=  2109  (  Pc  jUc  (uCx^s]  This  value  is  obvi°usly  greater 
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than  the  value  of  q  .  When  -fu,  increases  very  significantly, 
the  oxygen  on  the  surface  of  the  wall  is  not  enough  for  the 
total  oxidation  of  methane.  The  excess  methane  will  enter 
the  middle  section  of  the  boundary  layer  and  the  combustion 
within  the  boundary  layer  will  be  most  intense  and  there  will 
be  the  largest  amount  of  heat  released  from  this  process.  For 
example,  when  -f  = 0.533  the  local  maximum  temperature  within 
the  boundary  layer  is  ?300K,  exceeding  the  temperature  at  the 
outer  fringe  of  the  layer.  Different  from  the  phenomenon  of 
the  burning  of  methane  on  the  surface  of  the  wall,  this  combustion 
heat  will  not  be  transferred  directly  onto  the  surface  of  the 
wall.  The  heat  will  only  be  transferred  to  the  surface  of  the 
wall  after  a  certain  range  of  molecular  diffusion  and  heat 
convection.  By  increasing  the  amount  of  the  injection,  we  can 
thicken  the  boundary  layer  and  hence  reduce  the  heat  transfer 
onto  the  surface  of  the  wall.  By  combining  these  two  effects, 
we  can  see  that  the  heat  flow  actually  transferred  to  the 
wall  surface  will  be  slightly  lower  than  the  maximum  heat  flow 
generated  by  burning  on  the  surface  of  the  wall.  When  -fw  =0.533* 
the  heat  flow  can  be  43#  lower  than  the  qQ  associated  with  the 
case  of  no  injection.  Except  for  values  near  the  critical 
amount  of  injection,  the  heat  flow  decreases  roughly  linearly 
with  the  increase  in  the  amount  of  the  injection.  We  can  also 
say  that  the  decrease  in  heat  flow  is  proportional  to  (-fw)0*^®2. 
We  will  call  the  amount  of  the  injection  that  corresponds 


12 


\ 


to  the  maximum  heat  flow  associated  with  the  burning  process 
on  the  surface  of  the  wall  the  critical  injection  amount 
-fw*  .  Based  on  C-^(O)  s  C^(0)  =?  Cg(0)  =  0,  we  can  use  the 
algebraic  equations  and  the  boundary  conditions  to  solve  for 
the  values  of  the  other  concentrations  on  the  surface  of  the  wall. 

By  following  this  procedure,  we  have  obtained  C2(0)  and  -fw  ■=:  0.0556. 
The  result  is  shown  in  Figure  1  as  the  peak  of  the  solid  line. 


Figure  1.  The  variation  of  the  thermal  coefficient  on  the  surface 
of  the  wall  with  the  injection  parameter  -fw  in  the 
presence  and  in  the  absence  of  the  combustion  reaction 
associated  with  the  injection  of  methane  at  the  stagnation 
point  of  the  boundary  layer. 
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The  conditions  are  $=1,  Pr  =  0.71  and$(0). 

1.  methane  with  combustion  reaction  Le=1.4 

2.  methane  without  combustion  reaction  Le  *1.4 

3.  methane  with  combustion  reaction  Le  = 1 

4.  according  to  reference  £7]  and  an  estimated  Le=1.4 

5.  cal/kg 

2.  Assume  methane  is  an  inert  gas  for  the  basic  condition 
This  assumption  is  made  so  that  the  combustion  heat  release 

effect  can  be  examined.  We  set  C2  5  a  s  0,  eliminate 

equations  ( 8 ) — ( 10 )  and  then  solve  the  remaining  equations.  The 

dashed  line  in  Figure  1  is  a  result  of  the  monotonic  decrease 

of  heat  flow  with  -f  .  The  heat  flow  in  this  curve  is  low 

w 

comparing  to  that  in  the  solid  line.  The  difference  in  the  value 
of  the  heat  flow  between  the  case  in  which  combustion  exists  and  the 
case  with  no  combustion  reflects  the  effect  of  the  release 

of  the  combustion  heat.  The  difference  in  heat  flow  generally 

* 

increases  when  -f  increases  to  -f  .  Its  maximum  value  is 

w  w 

1090  rfcjUc(ucx)s]  °*^  ,  equivalent  to  60%  of  qQ.  As  shown  by 
the  solid  line  in  Figure  2,  the  difference  in  heat  flow  becomes 
smoother  with  any  further  increase  in  -fw.  The  dashed  line  in 
Figure  2  shows  the  ratio  between  heat  flow  associated  with 
combustion  and  that  associated  with  no  combustion. 

3.  Let  Le  ~1  in  the  basic  condition 
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The  last  term  of  the  energy  equation  vanishes  under 
this  assumption,  hence  making  the  equation  easier  to  solve. 

The  result  of  the  heat  flow  computed  for  this  situation  is 
shown  as  the  dot-dash  line  in  Figure  1.  When  -fw^0.1,  this 
heat  flow  is  close  to  the  solid  line.  Le = 1  has  very  small 
influence  on  the  equation,  the  type  of  temperature  cross 
section  is  similar,  the  difference  in  the  numerical  values 
is  small  and  the  ultimate  difference  in  heat  flow  is  not  high. 

All  these  indicate  that  the  special  case  where  Le  =1  can 
reflect  the  characteristics  of  the  cases  where  Le^l. 

4.  Let  the  temperature  on  the  wall  surface  vary  in 
the  basic  condition 

The  result  of  the  computation  is  shown  by  the  solid  curve  in 
Figure  3.  The  variation  in  heat  flow  is  less  than  6%  when  the 
wall  temperature  varies  from  1500K  to  3000K.  This  indicates 
that  the  influence  of  the  wall  temperature  on  the  heat  flow  is 
not  large. 


Figure,  2.  Differences  in  heat  flow  values  between  the  case 
Key:  1-  cal/kg 
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with  combustion  and  the  case  without  combustion,  with 
methane  as  the  injecting  medium.  This  figure  also  shows 
the  multiplication  factor  of  the  increased  heat  flow 
associated  with  the  situation  where  combustion  exists. 


Figure  3.  The  variation  of  the  heat  flow  on  the  wall  surface 
with  a  single  basic  parameter,  with  the  conditions 
that  the  basic  parameters  are  fixed  to  be  ^=1, 

Pr  =tO.?l,  Le  =1.4  and  0(0)  =  0.4and-fw  ♦  0.533 
1.  variation  with  0(0)  2.  variation  with  Pr 

3.  variation  with  Le  4.  variation  with  1? 

5.  Cal/kg 
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5.  Let  -f  *  0.533  in  the  basic  condition  and  observe 

W 

separately  the  effects  of  |  ,  Le  and  Pr 
Figure  3  shows  the  heat  flow  curve  associated  with 
0.5^1  <2  (double  dot-dash  line).  The  heat  flow  will  increase 
more  quickly  with  the  increase  in  C  .  We  can  say  that  it  is 
nearly  directly  proportional  to  (I.965  -  0.  744$  ).  This  result 
primarily  reflects  that  fact  that  the  increase  in  the  viscosity 
coefficient  will  cause  an  increase  in  the  frictional  coefficient, 
hence  affecting  the  energy  exchange.  The  dashed  line  in  Figure  3 
shows  the  heat  flow  associated  with  1  ^  Le  4.  The  heat  flow 

decreases  gradually  with  the  increase  in  Le,  nearly  directly 

0  2?4 

proportional  to  Le  *  .  This  result  mainly  reflects  the 

reduction  of  enthalpy  and  the  gradient  of  enthalpy  on  the  wall 
surface  while  the  temperature  gradient  increases  slightly.  The 
combined  effect  causes  the  heat  flow  to  decrease.  The  dot-dash 
line  in  Figure  3  shows  the  heat  flow  associated  with  0.5^Pr^l. 
The  heat  flow  decreases  rapidly  with  the  increase  in  the  Pr 
number,  nearly  directly  proportional  to  Pr0*®*4’.  This  result 
mainly  reflects  the  simultaneous  reduction  of  the  enthalpy 
and  temperature  gradients  on  the  surface  of  the  wall,  hence 
reducing  the  energy  exchange.  The  momentum  exchange,  however, 
increases  slightly. 

The  reason  why  we  are  considering  the  effect  of  the 
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individual  variations  of  2,  Le  and  Pr  on  the  heat  transfer 
is  to  try  to  compensate  for  the  drawback  that  the  real 
values  of  the  transport  coefficients  can  not  be  determined. 
After  obtaining  the  estimates  based  on  other  methods,  we  can 
estimate  the  heat  flow  by  interpolating  the  curves  shown  in 
Figure  3. 

The  discussion  of  the  conclusion  mentioned  above  remains 
to  be  examined  by  future  experiments  and  more  accurate  theories. 
The  transport  coefficients  also  need  more  further  studies. 
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